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Abstract 

The nonabelian screening potential is calculated in the temporal axial 
gauge. The Slavnov- Taylor identity is used to construct the three-gluon ver- 
tex function from the inverse gluon propagator. After solving the Schwinger 
- Dyson equation beyond leading order we find that the obtained momentum 
dependence of the gluon self-energy at high temperature does not correspond 
to an attractive QCD - Debye potential, but instead it is repulsive and power 
behaved (~ 1/r 6 ) at large distance. 
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1 Introduction 

In recent papers [1,2] Rebhan has investigated the screening (Debye) mass 
for QCD at finite temperature beyond the leading order, especially under the 
aspect of its gauge dependence (in the case of covariant gauges). This analysis 
has been performed in the framework of the Braaten - Pisarski resummation 
scheme [3] , and it is successful to define a gauge-independent screening length 
in next-to-leading order. However, it is logarithmically sensitive to the mo- 
mentum scale 0(g 2 T), e.g. to the nonperturbative gluon magnetic mass, and 
it differs from previously found results [4,5]. 

On the other hand possible modifications of Debye screening have been 
discussed in Refs.[6,7], where it is shown that the momentum dependence of 
the gluon self-energy may essentially affect the large distance behaviour of 
the potential. 

In this paper we calculate the static QCD coulor-singlet potential beyond 
the leading order by exploiting the infrared limit of the (time-time component 
of the) gluon self-energy n 44 (0, k). We work in the temporal axial gauge, and 
we use the full Schwinger-Dyson equation for n 44 , which includes the ring 
graphs as well as nonperturbative vertex corrections. Gauge covariance is 
guaranteed by the exact Slavnov- Taylor identities in order to construct the 
three-gluon vertex [8,9]. 

As a result we find strong indications for a power-like screening potential 
for high temperature QCD. 

2 The gluon polarization tensor and the three-gluon vertex 

The following calculations are performed in the temporal axial gauge [4], 
which is well suited to keep gauge covariance and to evaluate the necessary 
infrared limit. In this gauge, in which the gauge vector is parallel to the 
velocity vector of the medium, the thermal Green function technique is 
considerably simplified, especially since the exact gluon propagator has only 
non-vanishing spatial components, 



where we suppress the colour indices. X>y only depends on two scalar func- 
tions F(k) and G(k), which are related to the gluon polarization tensor U^ u 




(1) 
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as follows 



G(k) = \ (E n« + §n 44 ) , F(k) = ^n 44 . (2) 

A further important advantage of the axial gauge is provided by the exact 
Slavnov- Taylor identity for the three-gluon vertex function in the form 

r,rj c 7 (r, p, q) = ig f abc [D^(p) - V$(q)] , (3) 

where f abc denote the structure constants of the colour group. 

This equation being rather simple is used to construct the nonpertur- 
bative vertex function necessary for calculating the gluon-self energy in the 
infrared limit, namely the IT 44 (0, |q|) function, which is required to determine 
the static screening potential. 

The exact nonperturbative graph representation for this tensor is well- 
known (see e.g. Refs.[4,9]). In the axial gauge it consists of four diagrams 
for the pure gluon contribution. However, when one considers the IT 44 - 
component only two nonperturbative diagrams contribute (the more com- 
plicated ones are vanishing). At finite temperature T = 1/(3 the analytical 
expression for these two nonvanishing graphs is found to be 

n«(o,|q|) = lf^j^v M (4) 



•^"E / ^^( 2 P4)[Ai(p + g)rij 4 (p + g, -P, -q)T>ji(p)} , 



where Tf>l = -igf abc T ij4 . All funct ions Dij and Tjj 4 are understood to be 
exact ones. The main problem is to find the nonperturbative expression for 
the rY,- 4 (r, p, g)-function. 

Since we do not know the exact expression for the three-gluon vertex 
function, we try a nonperturbative ansatz, which we exploit in the following. 
First the vertex function is decomposed as 

r2£(?,r,p) = ri%(q,r,p) L + rt%(q,r,p) T , (5) 

and its transverse and longitudinal parts are constructed independently. In 
the following the longitudinal one is used in the form, which has been found 



3 



in Ref.[8]: 



r 4 t(g, r, p)^-^/ abc <!^(r4-p4) 



1 



/y 2 'p^' 



'G(r) 



I P 2 



p2 p2 



F(r) F(p) 



"P4 



[(pr)<Jij -Pirj-](r 4 -p 4 ) 
\ 1 F(r) 



q 2 — r 2 



F(r) r\ 



QiU(q-r) 3 



1 (F(p) F{q)\ 



(6) 



and which is assumed to be valid for any set of momenta including the soft 
ones. Eq.(6) satisfies the exact Slavnov- Taylor identity by using the standard 
inverse gluon propagator 

^( P )^(%-^)(p 2 + G( P )) + (l + ^)ri^, (7) 

and more details of its properties can be found in Refs.[9,10]. 
The transverse part of the ^^j(q,r,p)- function is not constrained by the 
Slavnov- Taylor identity of Eq.(3), but it is introduced to cancel the leading 
infrared singularities of the longitudinal vertex function (6). It is given by 



r£(?,r,p) T = -igf abc 



(pq)c^ +p j q l ( F(g) _ F(rY 
\ q 2 r 2 



q 2 — r 2 



r 4 



■ (rq^+r^ ( F(q) _ F(p)\ 



q 2 — p 2 



2 2 / 

q v ) j 



(8) 



and together with the transverse part of the function T^(q, r,p) 

nl c J (q,r,p) T = -tgf abc 
Eq. (8) is constrained by two identities 



-(pr)gj + (pqfo- ( F(q) F{ 
\ q 2 r 



qy 2 



rf w (g, r, p)r^ = , r, p)p^ = . 

4 



(9) 



(10) 



In constructing the vertex function T&j it is necessary to impose that 
it is well-behaved in the infrared limit, i.e. for g 4 = 0, q — > (and similar 
for the momenta p and r). Indeed, it is the case with respect to p and r 
for the longitudinal and transverse parts separately, but not for the infrared 
limit in q. However, the sum in Eq.(5) is infrared stable, as one can easily 
verify. Therefore imposing this condition not only motivates the inclusion 
of the transverse part Eq.(8), but it also fixes its magnitude and (minimal) 
structure. 

In order to determine n 4 4(0, |q|) both the longitudinal (6) and transverse 
(8) vertex functions are taken into account in Eq.(4). In the following we 
restrict the calculation of II44 only up to the 0(g 3 ) term. In this approxima- 
tion essential simplifications occur due to the different infrared behaviour of 
the functions G(p) and F(p) in the static limit. If required the function G(p) 
may be used as an infrared cutoff (and e.g. be identified with the magnetic 
mass), but in the final stage of the calculation we set G = 0. The dependence 
on the function F(p) is evaluated in the infrared limit by keeping only the 
p 4 = mode in the Matsubara sum of Eq.(4) after identifying the leading 
term of the n 4 4-function by n 44 = m 2 = g 2 NT 2 /3. The poles at p 4 = are 
regularized in the standard manner when using the temporal axial gauge by 
the prescription: J2 P4 l/p 4 = 0. 
After the lenghty algebra is performed we obtain 

TT fn Inh 9 ' N 4. 2 ' jlN I ^ I mV + ( P q)] 2 




2 + m 2 ] [(p + q) 2 + m 2 } 

2q 2 1 



q 2 p 2 J p 2 [(p + q) 2 + m 2 ] 



(11) 



where all integrals have no divergencies and can be explicitly evaluated. 
The final result to order g 3 becomes 

, , . v 2 g 2 NmT f 1 vr q 4 
n 44 (0, |q|) = m 2 + y — {-- + - — *— (12) 



+ 



2tt ( 2 8mW 
m 2 — q 2 (q 2 + m 2 ) 2 



m 



v^q 2 " 2m 3 



Vqg (q 2 + 2m 2 ) 2 V? 

arctan H — - _ arctan — — 

m 4m d vq 2 2m 



valid for real values of q. It is important to note that in the low momentum 
expansion of II 44 (0, |q|) a positive cubic term |q| is present in Eq.(12) - 
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instead of a negative linear one in the case of the one-loop approximation 
[11]. It originates from the product of the longitudinal - longitudinal modes 
of V tj in Eq.(4). 

At q = the static self-energy is 



n 44 (o) 



2 g 2 NmT 
m H 

47T 



(13) 



IE. + ±(fE)*/2 

3 4vr l 3 ' 



as obtained recently in Ref.[12]. 

For reference we quote here separately the 0(g 3 ) contribution due to the 
product of the longitudinal and transverse modes of V^, 



n^(0,|q|)=m 2 + 



g 2 NmT 



2ir 



m 



my 



-arctan[ 



m 



- 1 



(14) 



which coincides with the expression (in the covariant Landau gauge) derived 
by Rebhan [1,2] in the framework of the Braaten - Pisarski resummation 
scheme [3]. 

3 The screening potential beyond the leading term 

Applying linear response theory it is possible to relate 1144(0, |q|) to the 
gluon-electric potential. For the colour- singlet case it is represented by the 
Fourier transform [4]: 



V(v) = 



N 2 -l 



-9 



HT)J 



2N * y ~ ' J (27r) 3 q 2 + n 44 (0,|q| 
N 2 -lg 2 (T) 7 dqq sin(gr) 

+ n 44 (o,g)' 



(15) 



2N 



2ir 2 r J q- 




/ 



where r = |r| and q = |q|. 

Considering large distances r it is usually assumed that only the static 
1144(0, |q| = 0)-limit is essential for evaluating Eq.(15), and that in this limit 
the observable screening length is related to the electric mass m 2 E = 1144(0), 
as defined in the temporal axial gauge. Since the leading term for g — > of 
this quantity (named as the Debye mass) was found to be gauge-independent 
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it is the belief that the exponential Debye screening is also the result of 
multi-loop calculations, although the momentum dependence of the gluon 
self-energy is expected to have a rather complicated analytical structure, e.g. 
branch points. Even for the case that the only arising singularities are poles, 
it cannot be excluded that they are becoming complex ones as solutions of 

q 2 + n 44 (o, yq2) = 0, (16) 

such that oscillations in the r dependence of the potential are resulting. 

It was demonstrated in Refs.[6,7] that odd powers of y 7 *?; with the branch 
point at q 2 = 0, are responsible for the asymptotic power behaviour of V(r) 
for large distance r and for high temperature. 

In our case the situation is analogous to the one discussed in the paper 
[6], because of the term in IT 44 of Eq.(12), which is odd (cubic) in powers 
of v^l 2 "- 111 order to facilitate the comparison with the one-loop result the 
same notation as in Ref. [6] is used in the following. The convenient variables 
are q = mz , x = mr and t = g 2 NT/(8m) = 3m/(8T), and a dimensionless 
screening function S(x, t) is defined by 

N 2 - 1 a 2 (T) 



and 



*--'>-f/ i* + £/7£W (18) 



where the function f(z) is obtained from Eq.(12), 

„, . 2t r _ . 4 9 A . arctan(^) . 9 . 9 arctan(z/2), . 
f(z) = —[-1 + (1 - Az 2 - z 4 ) ^ + (z 2 + 2) 2 v 1 ' }. (19) 

7T Z IZ 

Applying contour integration (with the contour surrounding the first quad- 
rant in the complex z— plane) the leading term in the high (but finite) 
temperature limit becomes 

S(x t)~- 1 dyy*e-*y 
[ ' } -nJ [(l-y 2 + /(y)) 2 + tV/4]' [ } 

with f(y) = f(—iz) on the first Riemann sheet. Contributions arising from 
possible (complex) poles and from cuts with imaginary branch points may 
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be neglected, since they vanish exponentially. The integral (20) is evaluated 
asymptotically by keeping only the leading term for y ~ in the denominator, 
since the other contributions become exponentially small for large x, or are 
subleading for g — > 0. Thus a power-behaved repulsive potential 

V W | m _^^ ,2!) 

results. 

The dependence of the exact S(x, t) of Eq. (18) on x = mr for a fixed value 
of t, i.e. of the coupling g respectively, is plotted in Fig. 1 for t — 0.5 (Fig. 
la) and for t = 1.0 (Fig. lb). The dotted curve represents the exponential 
Debye form for t = at leading order, S(x, 0) = exp(— x). The full curve is 
the result of our analysis with f(z) of Eq.(19) and the dashed curve is from 
Eq.(14). 

One can see that after including the momentum dependence in H44 as de- 
scribed above the exponential behaviour of S(x, t) is significantly modified: 
screening is increased for values of x < 3 (that means a positive correction to 
the Debye mass), but at larger distances there is antiscreening due to the cu- 
bic momentum term in II44 with S(x,t) ~ — 24t/(nx 5 ). A similar behaviour 
has been discussed in Ref.[6] in the case of the one-loop approximation (for 
a comparison one may see Fig. 1 of Ref . [6] ) . 

The dashed curve in Fig. 1 corresponds to the case given by Eq.(14), 
which is analogous to the situation discussed in Refs.[l,2] (for the gauge 
parameter fixed by a = 0), but without introducing a magnetic term. No 
odd terms in the momentum q are present. The screening function S(x,t) 
tends to become oscillatory around its zero value ; this behaviour is more 
transparent for increasing values of the coupling as illustrated in Fig. lb: for 
this purpose a rather extreme value of t = 1.0 is chosen. For this large value 
even the full curve shows two oscillations, but for x > 7.0 it remains negative 
and power-behaved. The dashed one shows oscillations for all x > 10.0, but 
which are damped by exp(— x). 

In order to test the reliability of our calculation with respect to its infrared 
sensitivity we introduce a small magnetic mass term in the transverse part 



2 But it has been pointed out by A. K. Rebhan that after introducing a small magnetic 
mass term in Eq.(14) as described in [1,2] the corresponding function S(x,t) does not 
change its sign, i.e. it remains positive. 
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of the gluon propagator by assuming a non-vanishing, but constant value for 
the function G of Eq.(2). Numerically we find that the full curve in Fig. 1 
remains stable with respect to this modification. Indeed this should be the 
case, since the cubic term in the momentum dependence of n 44 does not 
depend on the transverse part of the gluon propagator, i.e. on the magnetic 
mass. 

4 Conclusion 

To summarize, we have analysed the static colour- singlet potential V(r) 
in QCD at high temperature. In the temporal axial gauge and beyond leading 
order V(r) turns out to be power-behaved (~ 1/r 6 ) and repulsive (Eq.(21)). 
Although the detailed shape of this potential is gauge dependent, the de- 
scribed qualitative features may actually be independent of the chosen gauge, 
i.e. as a result of the next-to-leading order treatment, the positive cubic term 
in q of 1144(0, q) already reflects the correct behaviour of the gluon self-energy 
for small momenta q ~ 0. We expect that this cubic term (contrary to the 
linear one appearing in the one-loop approximation) remains infrared stable 
and survives under the further resummations. So the stated problem turns 
out to be more complicated than for the usually considered ansatz with the 
Debye mass, and more detailed investigations are desirable as a future task. 
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Figure Caption 

Fig. 1 The screening function S(x, t) as a function of x for two values of 
t = 0.5 (Fig. la) and t — 1.0 (Fig. lb) as described in the text. The change 
in scales (logarithmic and linear) has to be noted. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9402361vl 



